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LONG-TERM  GOALS 

The  long-term  goal  of  the  project  is  to  meet  the  goal  of  ONR  DRI  NLIW,  which  is 
to  achieve  the  basic  science  understanding  that  leads  to  a  predictive  capability  that  will  be  able  to  tell 
when  and  where  non-linear  internal  waves  will  occur  and  what  effects  they  will  have  on  the 
hydrodynamic  and  acoustic  environment.  This  project  focuses  on  the  use  of  remotely  sensed  variables, 
together  with  models,  that  can  reproduce  and  predict  the  generation  and  structure  of  these  waves,  their 
evolution  during  propagation,  and  the  processes  controlling  dissipation. 

OBJECTIVES 

1)  To  determine  the  statistical  features  of  ocean  internal  waves  in  SCS.  Interpreting  ten  years  of 
satellite  synthetic  aperture  radar  (SAR)  images,  the  statistical  features  of  ocean  internal  waves  in  SCS 
will  be  determined.  The  statistical  items  will  include  the  wavelength  distribution,  distribution  of 
number  of  waves  in  a  wave  packet,  characteristic  half  width  distribution,  generation  location 
distribution,  occurrence  seasonal  distribution,  and  propagation  direction  spectrum  on  the  continental 
shelf.  The  statistical  analysis  includes  the  ocean  environment  conditions  and  its  seasonal  variability. 
The  results  will  provide  the  users  an  outline  of  internal  wave  behavior  in  SCS,  serve  as  a  basis  for 
empirical  prediction  of  internal  wave  behavior  in  SCS,  and  contribute  to  creation  of  a  predictive 
system. 

2)  To  understand  the  effects  of  topography/thermocline  on  the  evolution  of  solitary  internal  waves  in 
SCS.  Sharp  changes  in  bottom  topography  and  thermocline  depth  are  an  obvious  feature  of  the  SCS 
continental  shelves.  Combining  satellite  images,  field  observations,  and  solitary  wave  dynamics,  the 
effects  of  those  boundary  conditions  on  the  evolution  of  the  internal  solitons  will  be  examined.  The 
results  may  be  used  to  predict  the  behavior  of  internal  waves  propagating  over  the  shelf  slope. 

3)  To  explore  the  SAR  imaging  conversion  mechanisms  of  internal  waves.  On  the  basis  of  the  physics 
of  radar  imaging  of  ocean  processes,  the  mechanisms  for  polarity  conversion  (bright-dark)  of  internal 
waves  SAR  image  will  be  examined.  The  results  will  further  reveal  SAR  imaging  mechanisms  and  be 
used  for  SAR  image  interpretation. 
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APPROACH 


1)  Satellite  images  and  field  data.  Satellite  SAR  is  the  most  powerful  sensor  for  oeean  remote  sensing 
beeause  of  its  all  weather,  all  day  abilities,  and  high  spatial  resolution.  The  spatial  resolution  of  the 
state-of-the-art  satellite  SAR  images  reaehes  20  -  30  m,  and  the  swath  width  reaches  100  -  450  km. 
These  specifications  match  the  requirements  for  observing  the  ocean  internal  waves  with  a  spatial 
characteristic  scale  within  a  range  from  0.1  km  to  10  km.  The  SAR  images  used  for  the  project  were 
taken  by  the  ERS-1  and  2,  the  RADARSAT-1  and  2  satellites.  Historical  AXBT  data  will  be  used  to 
determine  the  vertical  thermal  structure  in  the  study  area.  The  data  will  be  obtained  from  NOAA 
archives,  which  are  in  the  public  domain. 

2)  Data  processing  methods.  Commercialized  image  processing  software  packets  will  be  used  for 
imagery  enhancement,  orthorectification,  filtering,  and  data  extraction. 

3)  Theoretical  analysis.  The  objectives  of  theoretical  analysis  are  to  understand  the  nature,  physics, 
mechanism,  and  laws  of  variation  of  the  studied  process,  to  derive  unknown  geophysical  parameters 
using  parameters  and  information  derived  from  SAR  image  interpretation  as  inputs,  to  analyze  the 
relation  between  the  studied  process  and  the  surrounding  environment,  and  to  predict  the  future 
development  of  the  studied  process. 

Dr.  Quanan  Zheng  serves  as  a  lead  PI  and  coordinates  the  project.  Zheng  will  focus  on  data  collection, 
satellite  image  interpretation,  and  the  theoretical  analysis.  Dr.  R.  Dwi  Susanto  from  LDEO,  Columbia 
University,  serves  as  a  Co-PI.  Susanto  focuses  on  image  processing,  statistical  analysis,  and  takes  part 
in  data  collection  and  the  theoretical  analysis. 

WORK  COMPLETED 

1)  Statistical  analysis.  The  statistical  analysis  of  IW  occurrence  in  the  northern  South  China  Sea 
(NSCS)  has  been  completed  using  seven  years  of  satellite  SAR  images  from  1995  to  2001 .  The 
statistical  analysis  of  sea  surface  boundary  conditions  has  been  completed  using  long-term  field 
measurements  of  sea  surface  wind  and  sea  state  in  NSCS  [Editorial  Group,  1982].  The  statistical 
analysis  of  IW  generation  boundary  conditions  in  the  Euzon  Strait  has  been  completed  using  20  years 
of  vertical  temperature  profiles  from  1980  to  1999  [Boyer  and  Eevitus,  1994]. 

2)  Dynamical  analysis.  Nonlinear  analysis  of  effects  of  shoaling  thermocline  on  the  IW  generation  has 
been  completed.  The  role  of  Kuroshio  playing  in  IW  generation  in  the  Euzon  Strait  has  been  analyzed 
using  linear  wave  models  and  the  Eourier  transform  methods.  Dynamical  analysis  of  bottom- 
topography-induced  stationary  IW  in  NSCS  has  been  completed. 

3)  Paper  preparation.  Based  on  the  results,  six  papers  (see  PUBLICATIONS)  have  been  prepared 
and  submitted  to  international-circulated  and  peel-reviewed  journals.  Three  of  them  are  accepted  and 
in  press. 

RESULTS 

1)  Statistics  of  IW  occurrence  in  NSCS.  Using  seven  years  of  satellite  SAR  images  from  1995  to 
2001,  statistics  of  IW  occurrence  in  NSCS  generate  the  following  results. 
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Latitudinal  distribution  of  IW  packets  inNSCS  is  shown  in  Figure  1  (partially  from  Zhao  et  ah, 
2004)..  Statistics  indicate  that  22%  of  IW  packets  distributed  in  the  east  of  1 18“E  obviously  originate 
from  the  Luzon  Strait,  and  78%  of  IW  paekets  west  of  I  I8°E  may  propagate  from  the  east  or  evolve 
from  the  solitons  originating  from  the  east  boundary  owing  to  the  fission  effect  of  shoaling 
thermocline. 


Longitude 


Figure  1.  A  map  of  latitudinal  distribution  of  IW  packets  in  NSCS.  Bold  lines 
represent  crest  lines  of  leading  waves  in  IW  packets  interpreted  from  SAR  images. 
The  rectangular  box  on  the  right  defines  an  IW  generation  source  region  for  the 

dynamical  analysis. 


Yearly  distribution  of  IW  occurrence  frequeneies  is  shown  in  Eigure  2.  One  ean  see  an  interannual 
variability,  the  frequeneies  in  1995,  1998,  and  2000  are  2  to  5  times  higher  than  that  in  other  years. 
This  interannual  variability  implies  that  there  are  long-term  and  large  seale  proeesses  modifying  the 
SAR-observed  IW  oecurrenee. 
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Figure  2.  Yearly  distribution  of  SAR-observed  IW  occurrence  frequencies  in  NSCS. 

Monthly  SAR-observed  IW  occurrence  frequencies  show  that  the  high  frequencies  are  distributed 
from  April  to  July  and  reach  a  peak  in  June  with  a  maximum  frequency  of  20%.  The  low  occurrence 
frequencies  are  distributed  in  winter  from  December  to  February  of  next  year  with  a  minimum 
frequency  of  1.5%  in  February  (Figure  3). 


Figure  3.  Monthly  distribution  of  SAR-observed  IW  occurrence  frequencies  in  NSCS. 


2)  Statistics  of  sea  surface  boundary  conditions.  Statistics  indicate  that  high  frequencies  of  low  wind 
conditions  (scales  0-1  and  2-3)  are  distributed  from  April  to  September  and  reach  a  peak  phase  from 
May  to  July  with  a  maximum  frequency  of  58%  in  July.  The  low  frequencies  are  distributed  from 
November  to  March  of  next  year.  The  frequency  distributions  of  high  wind  conditions  (scales  4-5,  6-7, 
and  8)  are  just  anti-phase  with  low  wind  conditions.  High  frequencies  are  distributed  from  October  to 
March  of  next  year  and  reach  a  peak  phase  from  November  to  January  with  a  maximum  frequency  of 
87%  in  December  (Figure  4). 

The  statistical  distribution  of  monthly  sea  state  frequencies  is  shown  in  Figure  5.  One  can  see  that  high 
frequencies  of  Wind  Wave  Scales  0-2  (for  the  wave  height  less  than  0.7  m)  are  distributed  from  April 
to  August  and  reach  a  peak  value  of  44%  in  August.  The  low  frequencies  are  distributed  from  October 
to  February  of  next  year  with  frequencies  lower  than  16%.  March  and  September  appear  as  transit 
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phases.  The  frequeneies  of  Seales  3-4  (for  the  wave  height  from  0.8  to  1.9  m)  have  an  almost  even 
distribution  with  a  mean  of  41 .8%  and  a  standard  deviation  of  2.5%.  The  distributions  of  Seales  5  (for 
the  wave  height  from  2.0  to  3.4  m)  and  6  (for  the  wave  height  from  3.5  to  6.0  m)  are  anti-phase  with 
Seales  0-2. 


Figure  4.  Statistical  distributions  of  monthly  mean  wind  frequencies  at  statistical  cells  centered  at 
22‘’N  IIS^E  in  NSCS,  The  numerals  marked  at  peak  bars  represent  wind  scale  ranges. 


Jan  Feb  Mar  Apr  May  Jun  Jul  Aug  Sep  Oct  Nov  Dec 


Month 


Figure  5.  Same  as  Figure  4,  but  for  monthly  mean  sea  state. 

3)  Nonlinear  analysis  of  effects  of  shoaling  thermocline  on  the  IW  generation.  The  study  proposes 
that  the  IW  generation  needs  the  neeessary  and  suffieient  eonditions:  initial  disturbanee  formation,  and 
wave  amplitude  growth.  Due  to  dissipation  effeet  on  the  disturbanee  energy,  only  fully  grown  waves 
have  a  ehanee  to  radiate  out  of  the  souree  region.  A  physieal  model  and  PKdV  equation  are  applied  to 
the  analysis  of  the  suffieient  eonditions  for  solitary  IW  amplitude  growth.  The  results  indieate  that  the 
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thermocline  shoaling  provides  the  foreing  to  soliton  amplitude  growth,  so  that  the  soliton  amplitude 
growth  ratio  (SAGR)  serves  as  a  deeisive  faetor  for  the  IW  oeeurrenee  frequeney.  Theoretieal  analysis 
prediets  a  linear  relation  between  the  two.  Applieation  of  theoretieal  models  to  field  measurements 
gives  a  eorrelation  eoeffieient  as  high  as  0.845  with  a  eonfidenee  level  of  99%  for  months  from  Mareh 
to  November.  The  linear  regression  gives  a  eorrelation  eoeffieient  (R2)  of  0.6519  and  a  SAGR 
threshold  (minimum)  value  of  0.90  for  IW  oeeurrenee.  Aeeording  to  the  theoretieal  solutions,  the 
eastward  propagating  disturbanees  have  no  ehanee  to  grow  up,  so  that  they  hardly  appear  in  the  eastern 
Luzon  Strait. 


4)  The  role  of  Kuroshio  playing  in  IW  generation  in  the  Luzon  Strait.  A  linear  model  deseribing 
western  boundary  eurrent  instability  generation  meehanism  for  the  oeean  internal  waves  is  developed. 
Analytieal  solutions  are  derived  from  a  zero-order  eomplex  frequeney  -  wave  number  relation.  The 
imaginary  part  gives  an  exponential  growth  rate  of  IW  in  the  form 


I  ^  in  whieh  k  is  the  wavenumber  of  IW,  g  is  the  gravitational 

aeeeleration,  f  is  the  Colorist  parameter,  and  V  is  the  meridional  eomponent  of  western  boundary 
eurrent.  One  ean  see  that  the  internal  waves  would  absorb  the  energy  from  the  western  boundary 
eurrent  and  then  grow.  The  real  part  of  the  relation  represents  a  dispersion  relation  of  the  generated 

internal  waves  in  the  form  of  ~  V^K/  ^  ^  ^  in  whieh  co  is  the  wavenumber  of 

IW.  A  signifieant  feature  of  this  dispersion  relation  is  that  the  wave  propagation  eharaeteristies  are 
elosely  assoeiated  with  the  instability  of  western  boundary  eurrent.  The  theoretieal  results  are  used  in 
the  ease  of  the  Kuroshio  east  of  the  Luzon  Strait.  Based  on  the  analysis,  it  is  found  that  for  the  western 
propagating  disturbanee,  the  Kuroshio  west  wing  is  unstable  and  the  east  wing  is  stable;  while  the 
reverse  is  true  for  the  eastern  propagating  disturbanee.  The  results  are  used  to  interpret  satellite  SAR 
images  of  the  oeean  internal  waves,  whieh  are  generated  in  the  Luzon  Strait  and  propagated  westward. 


=  /  dx)l  IdV  1 " 


5)  Dynamical  analysis  of  bottom-topography-induced  stationary  IW  in  NSCS.  The  satellite  SAR 
images  display  wave-like  patterns  of  the  oeean  bottom  topographie  features  at  the  south  outlet  of 
Taiwan  Strait  (TS).  Field  measurements  indieate  that  the  most  TS  water  body  is  vertieally  stratified.  In 
order  to  explore  the  meehanisms  and  to  determine  the  relations  between  the  SAR  imagery  and  the 
bottom  features,  a  two-dimensional,  three-layer  oeean  model  with  sinusoidal  bottom  topographie 
features  is  developed.  Analytieal  solutions  and  inferenees  of  the  momentum  equations  of  the  oeean 
model  lead  to  the  following  solutions.  1).  In  the  lower  layer,  the  topography-indueed  waves 
(topographie  waves)  exist  in  the  form  of  stationary  waves,  whieh  satisfy  a  lower  boundary  resonanee 
eondition  a  =  kCO,  here  a  is  an  angular  frequeney  of  the  stationary  waves,  k  is  a  wavenumber  of 
bottom  topographie  eorrugation,  and  CO  is  a  baekground  eurrent  speed.  2).  As  internal  waves,  the 
topographie  waves  may  propagate  vertieally  to  the  upper  layer  with  an  unehanged  wavenumber  k,  if  a 

frequeney  relation  satisfied,  here N2  and N3  are  the  Brunt-Waisala  frequeneies  of 

middle  layer  and  upper  layer,  respeetively.  3).  The  topographie  waves  are  extremely  amplified  if  an 
upper  layer  resonanee  eondition  is  satisfied.  The  theoretical  SAR  image  of  topographic  waves  is 
derived  on  the  basis  of  current-modulated  small  wave  spectra.  The  results  indicate  that  the  topographic 
waves  on  SAR  images  have  the  same  wavelength  of  bottom  topographic  corrugation,  and  the  imagery 
brightness  peaks  are  either  inphase  or  antiphase  with  respect  to  the  topographic  corrugation,  depending 
on  a  sign  of  a  coupling  factor.  These  theoretical  predictions  are  verified  by  field  observations.  The 
results  provide  a  physical  basis  for  quantitative  interpretation  of  SAR  images  of  bottom  topographic 
waves  in  the  stratified  ocean. 
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IMPACT/APPLICATIONS 


The  results  of  this  projeet  will  provide  the  users  a  statistieal  outline  of  internal  wave  behavior  and 
boundary  eonditions  in  SCS,  and  will  benefit  the  broader  oeeanographie  eommunity,  oeean 
engineering  industries,  underwater  navigation  and  operational  users.  The  results  may  also  serve  as  a 
basis  for  empirieal,  theoretieal,  and  numerieal  prediction  models  of  internal  wave  behavior  in  SCS,  and 
contribute  to  creation  of  a  predictive  system.  The  results  will  further  reveal  SAR  imaging  mechanisms 
and  be  used  for  SAR  image  interpretation. 

RELATED  PROJECTS 

The  PI  of  this  project  serves  as  a  CO-PI  for  an  ongoing  ONR  PO  project  titled  “Analysis  of  Fine 
Structures  of  Flows,  Hydrography,  and  Fronts  in  Taiwan  Strait”.  The  study  areas  of  two  projects  are 
immediately  adjacent.  Therefore,  two  projects  sometimes  share  the  same  data  resources  of  field 
observations. 
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